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ABSTRACT 


A reciprocating Gadolinium core in a regeneration fluid column in the 
warm bore of a superconducting solenoidal magnet is considered for magnetic 
refrigeration in 3.517 MW (1000 ton) applications. A procedure is 
presented to minimize the amount of superconducting cable needed in the 
magnet design. Estimated system capital costs for an ideal magnetic 
refrigerator of this type become comparable to conventional chillers as the 
frequency of reciprocation approaches 10 Hertz. A one-dimensional finite 
difference analysis of a regenerator cycling at 0.027 Hertz is presented 
which exhibits some of the features seen in the experiments of G.V. Brown. 
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NOMENCLATURE 


Column cross-section; Ampere 

Magnetic induction 

Magnetic induction at magnet center 

Magnet induction in central plane at wall of core 

Cost of magnet and magnetic metal core 

Cost of auxiliary equipment 

Specific cost of magnetic metal 

Specific heat at constant pressure 

Specific cost of superconducting cable 

Cost of magnet structural support 

Total cost of system 

Cost coefficient defined by eq. (16) 

Control Volume 

Form factor defined by eq. (21) 

Frequency of demagnetization (cycle frequency) 
Magnetic intensity 
Hertz (cycles/sec) 


Enthalpy 


h Specific enthalpy 

J Globally averaged current density in magnet winding 

J sa Characteristic current density of superconducting alloy 

K Unit of temperature, Kelvin 

k Thermal conductivity 

L Length of magnet winding 

M Magnetization 

M m Mass of core magnetic metal 

M s Mass of superconducting cable 

N Unit of force, Newton 

R.j Inner radius of magnet winding 

R 0 Outer radius of magnet winding 

P Cost term defined by Eq. (19) 

Pmin P at a min 

Ql Heat transfer to column end from source 

Qh Heat transfer from column end to sink 

0^ Refrigeration load rate 

6. Thermal input power 

in 

6 ^ Thermal output power 

out 


Specific refrigeration capacity of magnetic metal 
Specific entropy 

Temperature; Unit of magnetic field, Tesla 
Time 

Internal energy 

Speed of core relative to column 

Volume of magnet bore 

Speed of core relative to magnet 

Speed of column relative to magnet 

Speed of fluid in core, relative to column 
Volume of superconducting cable 

Volume of superconducting alloy in cable 

Volume of magnet windings including spacing 
Specific volume 
Unit of power, Watt 
Coordinate relative to magnet 


Coordinate relative to column 



Greek Symbols 


a Ratio of outer to inner radii for magnet winding 

/} Ratio of length to inner diameter for magnet winding 

finin Value of p at a m in 

8 Fraction of winding volume filled by superconducting alloy 

5 S Fraction of winding volume filled by superconducting cable 

5 sa Fraction of cable cross-section filled by superconducting alloy 
€ Porosity defined by Af/A 

~f Fraction of core filled by magnetic metal 

p Density 

p m Density of magnetic metal (non-porous form) 

p s Density of superconducting cable 

p 0 Magnetic permeability of vacuum 

Subscripts 
H high 

L low 

f fluid 

m magnetic metal 

min minimum 



I. INTRODUCTION 


Some materials become hotter (cooler) when magnetized (demagnetized). 
Refrigerators and heat pumps based on the effect can be imagined in principle. 
In practice, adiabatic demagnetization has been important in cooling below a 
few degrees Kelvin. This report considers the feasibility of commercial 
application of a magnetocaloric refrigerator operating near room temperature. 
The motivation is to find devices of greater capacity, economy or design 
flexibility than are found in conventional technology to meet particular 
applications. 

In 1976, Brown [1,2] at NASA/Lewis suggested the possibility of practical 
magnetocaloric devices at normal temperatures. The bulk availability of rare- 
earths such as gadolinium and the advent of higher field superconducting 
magnets with considerably less power consumption was intriguing because the 
magnetocaloric effect is stronger near the Curie point (Gd 293K) and with 
larger field changes. 

Brown discussed several possible thermodynamic cycles and selected the 
magnetic Stirling cycle with regeneration for further study. A proof of 
concept laboratory device was successfully demonstrated [3]. The commercial 
feasibility was thought to rest on economics and the thermodynamic performance 
of the regenerator. These are the factors addressed in this paper. 
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Comparisons are made here on the capital and operating costs of 
conventional 1000 ton chillers and a system using a reciprocating porous 
gadolinium core in a fluid regenerating column. 

Progress is reported on modelling the gadolinium - fluid column 
regenerator. This modelling is to eventually take into account rate dependent 
and irreversible processes, but hasn't done so yet. 

II. BACKGROUND 

Magnetic Cooling Devices 

The magnetocaloric effect was first observed by Weiss and Piccard [4] in 
1918. Prior to this Edison [5] and Tesla [6] had patented designs for 
refrigerators and engines based on the ferromagnetic-paramagnetic transition. 

Temperatures down to 1 K can be obtained to liquefy He. Debye [7] in 1926 
and Giauque [8] in 1927 independently suggested that lower temperatures could 
be produced by the adiabatic demagnetization of a paramagnetic substance. The 
method was successfully tested in 1933 by Giauque and MacDougall [9]. 
Adiabatic demagnetization has been used since then in low-temperature 
research. 
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In "one-shot" devices a paramagnetic substance is placed in thermal 
contact with a low temperature reservoir of He and a material to be studied at 
low temperature. An applied magnetic field causes thermal energy to flow from 
the paramagnetic material into the reservoir. Once thermal equilibrium is 
established the thermal contact to the reservoir is broken and the magnetic 
field is lowered to zero. The experimental material and the paramagnetic 
substance will then drop to a temperature below that of the reservoir. 

Refrigerators have been built using paramagnetic substances to maintain 
temperatures below 1 K for loads less than 1 mW [10, 11]. There is 
substantial interest in using similar devices to maintain low temperatures (1 
to 20 K) in superconducting devices and to cool instruments in space craft 
[13,' 14, 15]. These devices would need to handle loads greater than 1 W. 

There are basically two competing designs in current magnetic 
refrigeration research at temperatures above 1 K. One involves a porous 
magnetic material moving with a reciprocating motion in a fluid column. The 
other design uses a rotating wheel of magnetic material with a counterflowing 
fluid acting as the link between the source and sink. The two designs are 
referred to as the reciprocating and rotating designs, respectively. Both 
rotating and reciprocating magnetic heat pumps have been proposed for 
applications in space, laser amplifier cooling, helium 1 iquefication and 
industrial waste heat recovery [13, 14, 15, 16, 17]. 
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Reciprocating Magnetic Heat Pumps 


The reciprocating magnetic heat pump consists in part of the porous 
magnetic material, a fluid-filled regenerator column and the external magnet. 
A typical cycle is shown in Fig. 1. Here, the magnet is on continuously. 

In the mechanical cycle, process 1-2 includes magnetization of the 
magnetic material in the core. In the isofield process 2-3 the material is 
cooled. The cycle is completed by core demagnetization 3-4 and an isofield 
increase in temperature 4-1. In practice, the necessary translation would 
likely be vertical for proper fluid control. At steady operating conditions 
the fluid in the regenerator column is stratified with respect to temperature 
and has an overall temperature difference of T^ - T|_. In the column shown in 
Fig. 1 the left end is hot. The motion of the magnetic material with respect 
to the magnet is not always required. Instead, the magnet could be turned on 
and off with the magnetic material inside the bore. However, the nature of 
high-field superconducting magnets may favor the use of relative motion and a 
constant field. 

The details of the actual energy addition and rejection processes at the 
ends of the regenerator column are not shown. Various methods could be used, 
and these, along with the details of the relative motion and the spatial 
variation in the magnetic field, will determine what type thermodynamic cycle 
is achieved. Fig. 2 shows a cycle consisting of two isothermal and two 
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isofield paths which is the representation of the cycle considered here. If 
regeneration paths 2-3 and 4-1 were congruent as suggested by the dashed line, 
then the ideal cycle would have the same coefficient of performance as the 
Carnot device. (The dashed path could be achieved in practice by programming 
the partial magnetization of the core during warming, which is not considered 
initially. ) 

The reciprocating magnetic heat pump was first proposed by J. R. Van Geuns 
in 1966 [18]. More recently this device has been discussed in several papers 
by G. V. Brown [1, 2, 19]. Brown and Papell [3] have built and tested a small 
reciprocating magnetic device with adiabatic walls (no source or sink). The 
maximum field employed was 7 T producing a maximum temperature span of about 
80 K. In separate tests the lowest and highest temperatures attained were 241 
K and 328 K. Two factors which would limit the performance of an actual 
refrigerator were noted. The successful operation of a reciprocating 
refrigerator depends on the maintenance of a temperature gradient in the 
fluid, but the gradient in the test device was degraded by jets of fluid 
issuing from the core causing fluid mixing in the region behind it. Also, 
more surface area was needed to enhance heat transfer . 

Two other reciprocating magnetic refrigerators have been tested. Barclay 
et al. [20] built and tested a device which operated at source and sink 
temperatures of 2.2 K and 4.2 K, respectively. These limiting factors were 
noted: (1) frictional heating (mechanical contact), (2) viscous heating, and 
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(3) mixing owing to the motion of the porous core. C. Delpuech et al . [21] 
tested a double acting reciprocating magnetic refrigerator in 1981. This 
device has two paramagnetic cores and magnets. The cold section is located at 
the middle of the regenerator column and thermal energy is rejected from each 
end. The refrigerator was tested between 1.8 K and 4.2 K and produced nearly 
one-half watt of refrigeration capacity. 

Rotating Magnetic Heat Pumps 

The rotating magnetic refrigerator is arranged as a counterflow heat 
exchanger as illustrated in Fig. 3. The rim of the wheel is composed of 
porous magnetic material. Fluid is pumped through the porous rim as it 
rotates through high-field and low-field regions. 

A prototype to test the rotary magnetic heat pump principle was designed, 
built and tested in 1977 [22]. This device operated at room temperature and 
was a forerunner of a room temperature device that was reported on in 1981 
[23]. Also, a rotating magnetic refrigerator operating between about 2 K and 
4 K has been tested [24]. 

The two main problems with rotary designs are obtaining a high 
concentrated field at one location on the wheel and a zero field elsewhere, 
and controlling the flow of the fluid. These problems, along with lower than 
expected heat transfer between the fluid and magnetic material, resulted in a 
refrigeration capacity of 400 W, a Coefficient of Performance (COP) of 26 
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percent of the Carnot COP and a maximum AT of 7 K as compared to design goals 
of 1000 W, 70 percent and 40 K, respectively [23]. Barclay [25] has suggested 
that the flow problem might be alleviated somewhat by using a ferrofluid which 
would be driven through the porous material by magnetic forces. 

A Different Magnetic Heat Pump 

One other current device using magnetic materials, that does not fall into 
the previous two categories, is described in a 1984 patent by H. Nakagome and 
T. Hashimoto [26]. They envision a refrigerator composed of a magnetic 
material connected to one-way heat pipes. No published accounts of an 
operating device have been located. 

III. ECONOMIC ESTIMATES 


Capital Cost Model 

To obtain economic estimates the basic system of Fig. 1 with a steady 
magnet and a reciprocating magnetic porous core was chosen. Some simplifying 
assumptions were made. 

1) The porous core just fills the bore of the magnet windings. 

2) The vacuum field of the magnet is uniform in the bore with a value 
equal to that calculated for the magnet center. 
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3) Any metallic piece of magnetic core responds to the uniform applied 
field in the same manner as the center element of a long thin 
ellipsoid. Demagnetization is negligible. 

4) The magnet is assembled by stacking thin disk-shaped coils with 
adequate voids for cooling channels and structural members. 

5) The vacuum field at the magnet center is calculated using a global 
spatially averaged uniform current density. 

6) The refrigeration rate is proportional to the frequency of 
demagnetization of the core. 

7) Eddy currents are ignored. 

8) The coefficient of performance is the maximum COP: There are no 

irreversibilities. 

The logic that follows starts with a load specification. Then the amount 
of magnetic material to satisfy the load requirement is found. Finally after 
the superconducting cable is chosen, the shape of the magnet solenoid is 
optimized to provide the necessary field at minimum cable weight (for a 
uniform winding) . 

The refrigeration rate is 


\ \ V 


0) 


where q^ is the refrigeration capacity per unit mass of core magnetic metal 
per each demagnetization, 


M is the mass of the core magnetic metal, and 
m 
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f is the frequency of demagnetization. 


Therefore the required volume of the magnetic material is 

v . A_ 

m «L f ', 

where is the density of the magnetic metal in non-porous form, 

for the porosity of the core, 


( 2 ) 

Allowing 


where 


V b is the volume of the magnet bore, and 

7 is the filling fraction of the magnetic metal in the core. 


(3) 


The magnet is illustrated in more detail in Fig. 4 
where R 0 is the outer radius of the winding, 

Ri is the inner radius of the winding, 

L is the length of the winding, 

J is the globally averaged current density, and 
B 0 = B w = B is assumed. 

Defining 

R 

o 

a * R. 


it follows that 

V. = 2?r R 3 0. 
b i 


(4) 

(5) 

( 6 ) 
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The windings and their spacing occupy a volume, V w given by 

V w « (a 2 - 1) V b . (7) 

The superconducting cable has a volume V s 

V s = 5 S V w (8) 

where S s is the fraction of the winding volume filled with superconducting 
cable. 

A typical superconductor cable consists of tiny filaments of 
superconducting alloy embedded in copper. For the purposes here, it is 
necessary to specify the cable in more detail. Let 6 sa be the fraction of the 
cable cross-section which is superconducting alloy, then 

Vsa « *sa V s (9) 

where V sa is the volume of the superconducting alloy. 

Defining S S S sa 6 S (10) 

then V sa = SV W . (11) 

Now the global average current density J can be related to the known 
superconducting alloy characteristic current density, J sa . 

J =5 d s a (12) 

A complication arises in that q. and J each depend on the field. 

L sa 


- J (B) 
sa sa 

(13) 

q L = q L (B ’ V 

(14) 
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Thus one must choose B, then find the volume of magnetic material, V m , through 
equation (2), then proceed to find the optimum winding shape starting with 
equations (13) and (12). 


Since V m and V sa (expected to be high cost variables) are coupled by B one 
may seek to estimate the system cost as a function of B. 

Let Ct be the total cost of the system 

C t * C m M m + C S M S + C s t + C a (15) 

where C m is the specific cost of the magnetic metal, 

is the specific cost of the superconducting cable, 
is the cost of the magnet structure, and 
is the cost of the auxiliary equipment. 


Cm 

Cs 

Cst 

Ca 


Assume 

Define 

where 
form, and 


Pm 


C s t + C S M S - Ci C S M S . (16) 

C s C m p m V m + CiC s /> s V s (17) 

is the density of the magnetic material in non-porous 


p s is the density of the superconducting cable. 

Rewrite equation (17) using equations (3), (7) and (8): 

C * Cm/>[nV m [1 + P] (16) 


where 


P = 


'C.C p & 

1 s s s 


C pi 
mm 


(« - 1 ). 


(19) 
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Cost Minimization Logic 


Recall equation (2). 

v -A. 

m f Vm 

The application fixes Q|_, Tl. 

Choose the frequency, f, and the magnetic metal, p m . 

Choose B which fixes q|_ = qL (B, Tl). 

Now V m is fixed in equations (2) and (18). 

Choose the cable: C s , p s , 6 sa . 

Choose the filling fractions: y, 8 S . 

B also fixes J$ a = dg a (B). 

Now one can minimize the cost C of the magnet and magnetic material for given 
field B by minimizing Cj (a 2 - 1) in equations (19) and (18). 

Assuming that Cj is not a sensitive function of the design variables, one 
seeks to minimize a for a given B. 

Minimum Radius Ratio, a 

Following Reference [28] for this type magnet we have 

B = J Ri F (20) 
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where F is a form factor 


F = F 



f 2 2 

ct + (ft + ) 

1 + (1 + /) 


1/2 1 

172 “ 


( 21 ) 


where /x 0 is the magnetic permeability of vacuum. 

In equations (20) and (21), B and J have been fixed. The parameters ft, /J 
relate only to the shape of the magnet: R 0 , Ri, L. One seeks to find the 
shape which minimizes ft subject to the constraints of equations (20) and (21). 
The latter can be recast in an informative manner. Combine equations (3), 
(6), (12), (20) and (21) to get: 


31 

6 fv 

, J 
| sa 

UttTJ [ 

mj L B J 


or, using equation (2) 


s 

M 


[*l1 


1 1 


J 1 
sa 

F 3 ' 

[2tt7 


If 


q LV 


B 


( 22 ) 


(23) 


Since the parameters on the right have been fixed, 

D 

Hr = constant. 

i- 


(24) 


From equation (21), 0/F 3 is a function only of ft and 

One can easily find ft m i n subject to the constraint in equation (24). A 
typical curve for ft vs /} is shown in Fig. 5. Now let j9 m -j n stand for the value 
of /} associated with a m in- 
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Optimum Magnet Dimensions 


Combine equations (3) and (6) after inserting the values for a m i n , p m \ n to 


get 


m 

2 n -y /J . 
min 


1/3 


From equation (4) 


R = a . R , . 
o min i 


The "build" of the winding defined as R 0 - R^, is 

R o ' R i “ R i ( a min ' 1) 

From equation (5) , 

L = 2 0 m i n Ri . 

Of course the inner and outer diameters are 

Di = 2 Ri 
^o * 2 Rq • 


Input Data Case I 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 


For the parameters of equation (23), 

Load Requirements: 

Ql - 3517 kw (1000 Tons of Refrigeration) 
f = 1Hz 
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Geometric Design: 


7 * .8 

5 « .1125 (<5 S = .9, 6 sa = .125) 

Magnetic Metal: Gadolinium 

p m = 7.9 x 10 3 Kg/m 3 

C m = $200/Kg 

q L = [ . 589B - .0817 B 3y/2 ] KJ/Kg at T L = 280K 

B in Tesl a, Ref. [29] , 

Superconducting Cable: 

1 part superconducting alloy, 7 parts copper 
p s = 8.6 x 10 3 Kg/m 3 

C s = $66/Kg, Ref. [30]. 

Superconducting Alloy: 

J sa = J S a( B ) as given below. 

J sa = (5.90) 10 10 (.773 ) b Amp/m 2 , 

for B in the range 2.5T < B < 12T. Ref [31]. 


Cost Formula - Equations (18) and (19): 
C m Am = 1-58 x 10 6 $/m 3 
C s p s - .568 x 10 6 $/m 3 . 



Equation (18) becomes 


C = $ ( 1 . 58E6 ) [1 + P min ] V M , 


(31) 


Where P m i n is given by setting a = a m i n in equation (19). 


Equation (19) becomes 


P = ( . 404 ) C i (a m in 2 - !)• 
In the above, the cost of the Gd is: 

$(1.58E6)V M 

The cost of the superconducting cable is: 


$ ( 1 . 58E6 ) V M 4 • 
L 1 


The cost of the superconducting magnet is: 

$(1.58E6)V M P. 


(32) 


(33) 


(34) 


(35) 


Input Data Case II 


The only change made for the second evaluation was to increase the 
frequency by a factor of 10. 
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Results 


The results of the calculations for the costs of the Gadolinium and the 
cable are given in Tables 1 and 2, along with estimates for the costs of the 
magnet and the auxiliary equipment [30]. 

The magnet cost estimates have a large uncertainty because of the lack of 
appropriate data. There are commercial warm bore magnets for the Magnetic 
Resonance Imaging and other applications of about the same bore diameter and 
length as those in Tables 1 and 2. The MRI costs are strongly escalated by 
the need for field homogeneity over a relatively large volume. Field 
homogeneity is probably not so critical in magnetic heat pumps. On the other 
hand, the cost of the latter will escalate due to the need of the magnet to 
withstand large internal forces as the Gadolinium is withdrawn from the 
magnet. 

For the cost of the magnet, C s + C s t> we have simply used the expression: 
C s + C t = $(14B + 19) k , B in Tesla. 

This cost formula fits fairly well to marketed non-MRI magnets with bore 
shape and size similar to those of Tables 1 and 2. The $ 1 9k represents the 
costs of a rack of electronics and a closed-loop liquid helium refrigeration 
system, both used for the operation of the magnet. The costs reflect small 
quantity production and do not include the significant reductions expected 
from high quantity manufacturing. 
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Table 1. Capital Cost Estimates Vs. Applied Field, f = 1 Hz. 


3T 


Vm> m 3 


.3315 

j3/F 3 , A 6 /N 3 


7.036E25 

®min 


1.00265 

0min 


1.4 

Ri , cm 


36.12 

R 0 -R-j , cm 


.0957 

L, cm 


101.1 

Cost, Gd $ 


524k 

Cost, Cable, $ 


1.1k 

Cost of Magnet, 

s 

61k 

Cost of Auxiliaries 

30k 

Cost of System, 

$ 

615k 


61 

21 

121 

.1908 

.1438 

.1212 

4.970E23 

1.093E22 

3.819E20 

1.01384 

1.04972 

1.15458 

1.4 

1.45 

1.5 

30.04 

27.02 

25.24 

.416 

1.112 

3.90 

84.1 

78.4 

75.7 

301k 

227k 

192k 

3.4k 

9.4k 

25.8k 

103k 

145k 

187k 

30k 

30k 

30k 

434k 

402k 

409k 


Table 2. Capital Cost Estimates vs. Applied Field, f - 10 Hz 


3T 


V m , m3 .03315 

)3/F 3 , A 6 /N 3 7.036E24 

°r m i n 1.00572 

0mi n 

R-j, cm 16.76 

R 0 -R-j, cm .0959 

L, cm 46.9 

Cost, Gd, $ 52.4k 

Cost, Cable, $ .2k 

Cost of Magnet, $ 61k 

Cost of Auxiliaries, $ 30k 

Cost of System, $ 143k 


61 

21 

121 

.01908 

.01438 

.01212 

4.970E22 

1.093E21 

3.819E19 

1.02991 

1.10808 

1.3425 

1.45 

1.5 

1.65 

13.77 

12.39 

11.34 

.419* 

1.339 

3.89 

39.9 

37.2 

37.4 

30.1k 

22.7k 

19.2k 

.7k 

2.1k 

6.2k 

103k 

145k 

187k 

30k 

30k 

30k 

163k 

198k 

236k 
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Separately listed as auxiliary equipment is the mechanical system required 
to drive the motion of the regenerator and the magnetic core. If the ideal 

magnetic refrigerator operates between 280K and 310K, the ideal COP is 9.33, 

and the work required for a 35.17 MW load at f = 1 Hz would be 377 kJ of work 
per cycle. If the stroke length were 1 m, then the average force would be on 
the order of 377 kN with a peak force about twice as high. At f = 10 Hz, the 
average level would be about 37.7 kN. Since the actual motion would be 
vertical, the weight of the Gadolinium must be added, which increases the 
estimated level about 8%. 

Even the case with the greatest forces (B = 3T, f = 1) can be met with 
hydraulic piston/cylinders and electromagnetic actuators. A non-magnetic 
stainless steel rod of say 5 cm diameter would suffice as the major driver 

(one rod on each end of the metal core). If the rods had to traverse the 

fluid column, each would use only a small percentage of the volume available 
(an effect neglected in Tables 1 and 2). However, by using some kind of 
internal latching between the core and the regenerator walls in paths 1-2 and 
3-4 (Fig. 1), it may be possible to have the major drivers attached to the 
ends of the regenerator column. 

Another set of piston/cylinders, valves and rods (the minor drivers) are 
needed to move the column relative to the core in paths 2-3 and 4-1. The 
mechanical power requirements here are much smaller than for the major 
drivers. The initial costs estimated for the mechanical system are [39]: 
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$ 6k 
4k 
6k 
4k 
10k 

TOTAL $ 30k 


two cylinders 
two pumps 
valves 
controller 
miscellaneous 


Operating Cost 

An ideal thermodynamic cycle has been assumed. The quantitative effects 
of irreversibilities due to core/fluid interactions, eddy currents in the 
Gadolinium and other causes have not been assessed yet. The irreversibil ites 
do increase with cycle frequency and will offset to some extent the capital 
cost advantage of higher frequency observed in Tables 1 and 2. 

A serious loss occurs in the mechanical system in the conversion of motor 
shaft power to cylinder rod power. In conventional systems this conversion 
efficiency is about 80% [39]. However, the existing trend in hydraulics to 
much better efficiencies by using higher pressures and very efficient pumps, 
motors, and valves looks very encouraging [39]. 
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rV. REGENERATOR COLUMN 


The conceptual design selected for detailed system analysis is the 
reciprocating core in a regenerative fluid column within the bore of a steady 
superconducting magnet as in Fig. 1. The analysis is the subject, in part, of 
a dissertation in progress by one of the authors [32]. The first task was to 
treat the regenerator fluid and core using a simplified one-dimensional, 
transient model. 

The assumptions in the model are: 

1. The core is assumed to have a porous structure composed of gadolinium 
having a uniform porosity. 

2. Temperature gradients in the fluid and core normal to the direction 
of motion are negligible. The fluid and gadolinium in the core are 
in thermal equilibrium in any cross-section. 

3. Viscous forces and inertial forces are ignored. 

4. Fluid properties are independent of temperature. 

5. The magnetic field intensity, H, is a known function of position and 
steady in time, and is taken to be the vacuum field of the magnet. 

6. The entropy of the gadolinium is a known function of temperature and 
magnetic field intensity. 

7. The gadolinium is rigid and the fluid is incompressible. 
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8. The velocity of the core with respect to the regenerator column is 
constant during core traversals in the column. 

9. There is no dwell time. Either the core or column (or both) is in 
motion relative to the magnet at all times. 

10. Any effects of magnetically induced eddy currents are negligible. 

The thermodynamic properties of gadolinium are given by Griffel [33], 

Brown [2], and Benford and Brown [34], The general thermodynamic relations 
for magnetic materials are given by Hatsopoulos and Keenan [35] and Booker 
[36]. The relation between the applied fields (no magnetic material present) 
and the internal fields in the gadolinium in place is taken to be that of an 
ellipsoid of gadolinium [37,38] with no demagnetization effect. 

To write the energy balance, fix a reference frame (x) to the left end of 
the column in Fig. 1 and assume the Gd core is moving to the right with speed 
V. For a differential control volume located at x in this frame 
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Net Enthalpy Input Rate = 
where U is the internal energy, 
t is the time, 
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0 is the thermal power, 

W is the rate of magnetic work, 

H* is the enthalpy, 

h is the specific enthalpy, 

m, f are subscripts for magnetic metal, fluid, 

p is the density, 

A m? f is the cross-section area of metal or fluid, 
k is the thermal conductivity, 

T is the temperature, 

H 0 is the magnetic permeability of free space, 

H is the magnetic intensity, 

M is the magnetization, 

V is the metal speed, relative to the column, and 

Vf is the fluid speed in the core, relative to the column. 


f V f (40) 
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In this model the two magnetic vectors are col inear. 

Combining the equations, dividing by AAx, and taking the limit, one 
obtains 
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where A is the column cross-section area, 

v m is the specific volume of the solid, 

e is the porosity defined by € s Af/A, 

Cp is the constant pressure specific heat, 

and use has been made of the relations: 

dh = C dT 
P 

6 —-- d - ,M d - 
dt at 3x 
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It is convenient to change variables because the thermodynamic properties 
of Gd are available [33,34] in the form of the entropy function 

s = s(T,P,H) 

at atmospheric pressure. From Ref. [35], 

Tds = dh - /x 0 vHdM -vdP (43) 

so that for the metal in a constant pressure process, 
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which is the factor to be transformed. 
From calculus, at constant pressure, 
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By this transformation of variables, the energy equation becomes, after letting 
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is the total change in the field intensity observed at the magnetic material. 
The first term appears because the field appears to be time varying relative to 
the column due to the column motion. This term can be referenced to the magnet 
fixed coordinate (X) which removes the time varying component. Then 


d(vH) 3(vH) 
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ax 
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where V L = V , + V. V . is the speed of the column relative to the magnet, 
abs col col 


If e = 1, there is no Gd and the equation is 
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as expected, then the energy change is due only 
The boundary conditions are: 

1. At the regenerator column ends: 
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to conduction in the fluid. 


Initially adibatic, — = 0; then later isothermal, T = fixed. 
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These conditions represent the continuity of temperature and heat flux. 
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Initially, the core and regenerator column are in thermal equilibrium at 
the temperature of the environment. The column ends are adiabatic so that the 
mechanical cycle of Fig. 1 will cause a temperature gradient to develop in the 
regenerator column, the left end being the hotter. The thermodynamic path is 
not a cycle because after each mechanical cycle, the temperature profile 
differs. However, eventually a thermodynamic cycle should result as the axial 
thermal conduction in the fluid limits the maximum temperature difference 
between the column ends. 

After the column ends have reached temperatures suitable for refrigeration, 
the column ends are to be put into appropriate thermal contact with the source 
and sink of the refrigeration scheme. As a step in that direction, the first 
calculations included thermal reservoirs in contact with perfectly diathermal 
walls at the column ends. Heat transfer with the reservoirs occurs via 
conduction in the fluid and Gd. 

Solution of the energy equation was implemented on a computer using a 
finite difference scheme. Non-dimensional izing the equation was attempted. 
This proved to be very difficult due to non-constant coefficients and the lack 
of global geometric scales that can be used for references. Therefore a 
dimensional approach was used initially. 

Some implementation difficulties, discussed more completely later [32], are 
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categorized as: 

(1) Node Types. 

The finite difference nodes are picked to be fixed relative to the 
regenerator column. As the core moves along the column the character of 
the nodes changes. Also, the core-fluid boundary is in general not at a 
node location. The numerical scheme must recognize different node types 
and use various schemes to calculate new temperatures. This leads to 
bookkeeping problems with the nodes. 

(2) Stability and Convergence. 

This is the usual problem with the ratio of step sizes. The space 
increment is chosen arbitrarily and the time increment At is computed to 
keep the coefficients in the numerical scheme positive. The increment At 
is also checked so that the core advances through the column by a set 
fraction of a space step. 

(3) Convective Terms. 

The convective terms jthose with ^ j had to be replaced with a one-sided 

upwind difference to improve stability. As V, the relative velocity between 

3T 

core and column changes sign, the — terms change relative to the upwind direc- 

dX 

tions. A higher order difference was also tried but did not improve stability. 
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Results 


A. Adiabatic Column Ends, Initial Thermal Equilibrium 

Some initial predictions from the model are given in Figs. 6-8. The 
component sizing was chosen independently of that in Section III. The 
regenerator column has a length of lm, the core, 0.2m. The fluid is half 
water, half methyl alcohol by volume. The porosity of the core (open 
area/section area) is 0.80. The cylinder and ends of the regenerator are 
adiabatic. The vacuum field of the superconducting magnet is assumed constant 
in time and given by American Magnetics (see Appendix) for their 8 Tesla unit 
operating at 6T maximum. 

At the start of the first cycle, the fluid and gadolinium are uniform in 
temperature at 295K. The end of the regenerator nearest to the magnet is lm 
from magnet center and the core is near that end. After a half-cycle, the core 
and magnet are concentric, and the core is near the other end of the 
regenerator. The speed of the column relative to the magnet is O.lm/s or zero. 
Also, the speed of the core relative to the magnet is O.lm/s or zero. The cycle 
period is 37 seconds with no pauses and with velocities given by appropriate 
step functions. 

Fig. 6 shows the temperature profile within the regenerator after 14.5 and 
15 cycles. During this computer run, the minimum separation between the core 
and regenerator end was 5cm at each end. The column position, 0.00, marks the 
near end of the column at the start of a cycle. The drop in temperature at the 
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other end between the two profiles shows the cooling effect of removing the 
gadolinium from the magnet (followed by repositioning the column to complete 
the mechanical cycle). At present the model includes thermal conduction in the 
column, but excludes mixing. The gradients at the column ends reveal that 
axial conduction has a minor effect. 

Fig. 7 shows the maximum and minimum temperatures in the column at each 
half-cycle. The nearly isothermal sections again reveal the minor effect of 
axial conduction. The curvatures suggest that asymptotes will be approached, 
but Th-Tl is already close to the value found in the experiment by Brown and 
Pappel [3]. 

Fig. 8 gives the results of a run in which there is no dead space between 
the gadolinium and the ends of the regenerator. The temperatures at the column 
ends are now changing by a larger amount with each cycle. 

B. Isothermal Column Ends, Initial Internal Thermal Gradients 

The column is operated adiabatically as before until the highest 
temperature in the column goes above 310K and the lowest temperature in the 
column goes below 280K. Then reservoirs at 310K and 280K are coupled 
di athermal ly to the hot and cold ends, respectively. This is done 
computationally by holding the end nodes of the column at 310K and 280K. The 
core is assumed to travel the entire column length so at various times the core 
transfers thermal energy directly to the reservoirs. The heat transfer is 
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calculated from the temperature gradient and conductivity at each end for each 
increment of time. This is summed over a mechanical cycle to get the energy 
into the cold end and the energy rejected at the hot end. Once a thermodynamic 
cycle is achieved, the difference in these two quantities equals the work into 
the system. 

Figure 9 shows the energy exchange with the reservoirs per mechanical 
cycle versus the number of successive mechanical cycles. The porosity is now 
0.50. Qh is the heat transfer at the sink per mechanical cycle, and Q|_ is at 
the source. It appears that asymptotic values may occur, thus signalling the 
achievement of thermodynamic cycling. 

In Fig. 10, the values of 



are plotted versus the number of successive mechanical cycles. The denominator 
does not represent the work involved because a thermodynamic cycle has not yet 
been obtained. However the data of Figs. 9 and 10 are not inconsistent with an 
asymptotic approach to a COP in the vicinity of the ideal COP of 9.33 for a 
Carnot refrigerator with reservoirs at 310K and 280K. 

Such a limit will not be reached in this model because different sections 
of the Gd execute different thermodynamic paths. The thermodynamic paths are 
shown in Fig. 11 for three sections (left end, middle, right end) for the 30th 
mechanical cycle. 
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The program is stopped arbitrarily when both Qh and Q|_ in the Nth cycle 
differ by less than 1% from their respective values in the (N-l)th cycle. The 
same criteria for program cut-off was used in calculations starting with 
different core porosities. The results are given in Fig. 12, where 


Q 


L 


Q H 


- Q 


L 


is plotted versus core porosity. Any conclusions should be drawn cautiously 
since Figs. 9 and 10 suggest a closer approach to the asymptote may be needed. 
The result at low porosity may not be in error even though it exceeds the 
Carnot limit of 9.33, because a thermodynamic cycle has not been achieved. 
Clarification of this was not pursued because research has been started on a 
two-temperature model . 


V. DISCUSSION/CONCLUSION 

Large conventional chillers (1000 ton = 3.517 MW) have a capital cost near 
$ 1 50k , an electric operating use near 600 kW, an actual Coefficient of 
Performance of 5.86 and an ideal COP of 10 when operating between 280K and 
308K. 

The estimated system capital costs for ideal magnetic refrigeration 
systems of the same capacity become comparable as the cycle frequency 
approaches 10Hz, as shown in Table 2. 
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A significant counter trend is that the irreversibilities in the 
regenerator due to the relative motion of the fluid and the porous metal 
increase with cycle frequency. Such irreversibil ites have not yet been treated 
in our modelling of the regenerator. 

The capital costs for the magnetic refrigeration system appear to be 
driven by the cost and temperature-entropy characteristics of the magnetic 
material (Gd in this report) and by the structure/assembly costs of the magnet. 
Neither have been optimized. Magnetic materials of slightly less performance 
and considerably less cost are available and may be suitable. No attempt has 
been made to incorporate cost reductions which accompany high quantity 
manufacturing. 

Future effort toward reducing the uncertainties in the magnet cost 
estimate may not be warranted until estimates of the COP of actual regenerators 
have been completed. The appropriate literature to do so seems scarce. In any 
case, an experimental program would be necessary. 
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Column Temperature Profile at 14.5 and 15 Cycles 
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Maximum and Minimum Temperatures in the Column versus Time. 
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Column Temperature Profile after 14.5 Cycles 
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Heat Transfer Per Cycle for Successive Mechanical Cycles 
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T vs. S of Gd Core for Mechanical Cycle No. 30 




APPENDIX 



The magnetic field profile for the modelling of the regenerator is given 
below. The profile for an 8T magnet supplied by American Magnetics was 
arbitrarily multiplied by 3/4. The values are for the axis locations. 

Distance from center, m Field, Tesla 

0.0000 6.000 

0.0127 5.966 


0.0254 


5.850 


0.0381 

0.0508 

0.0635 

0.0762 

0.0889 

0.1016 

0.1143 

0.1270 


5.700 

5.275 

4.669 

3.776 

2.729 

1.829 

1.210 

0.818 


0.1397 

0.1524 

0.1651 

0.1778 

0.1905 

0.2032 

0.2159 

0.2286 

0.2413 

0.2540 

0.2667 

0.2794 

0.2921 

0.3048 


0.572 

0.413 

0.308 

0.236 

0.184 

0.147 

0.119 

0.098 

0.082 

0.0691 

0.059 

0.051 

0.045 

0.038 
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